A model for net-baryon rapidity distribution 
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In nuclear collisions, a sizable fraction of the available energy is carried away by baryons. As the 
baryon number is conserved, the net-baryon B — B retains information on the energy-momentum 
carried by the incoming nuclei. A simple and consistent model for net-baryon production in high 
energy proton-proton and nucleus-nucleus collisions is presented. The basic ingredients of the model 
are valence string formation based on standard PDFs with QCD evolution and string fragmentation 
via the Schwinger mechanism. The results of the model are presented and compared with data at 
different centre-of-mass energies and centralities, as well as with existing models. These results show 
that a good description of the main features of net-baryon data is possible in the framework of a 
simplistic model, with the advantage of making the fundamental production mechanisms manifest. 
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I. INTRODUCTION 

In hadron-hadron, hadron-nucleus and nucleus-nucleus 
interactions a sizable fraction of the available energy in 
a collision is carried away by baryons [l], 0, 9 ■ As the 
baryon number is conserved, the measured net-baryon, 
B — B, keeps track of the energy- momentum carried by 
the incoming particles. An important question to be 
asked is: how does the fraction of energy carried by the 
net-baryon evolve as a function of the centre-of-mass col- 
lisional energy per nucleon, y/s ? This question is im- 
portant because, on one hand, a decrease of the fraction 
of energy going into the net-baryon implies more energy 
available to create the deconfmed quark-gluon state of 
matter, and on the other hand, such a decrease may re- 
duce the possibility of producing fast particles in very 
high energy cosmic ray experiments. For more than 30 
years, since the ISR at CERN, particle production stud- 
ies have been limited to mid rapidity. Fortunately, with 
RHIC, large rapidity data became available, and, hope- 
fully, the same will happen with the LHC. In fact, if one 
does not measure the physics at high rapidity the most 
elementary physical constraint, namely energy conserva- 
tion, cannot be applied 0, Q . 

In most of the existing Monte Carlo models @, 0, H, , 
the physics of net-baryon production is very much ob- 
scured by the complexity of extensive and detailed codes. 
Often, the basic production mechanisms do not appear 
in a transparent way. In this paper we present a simple 
and consistent model for net-baryon production in high 
energy proton-proton (p-p) and nucleus-nucleus (A-A) 
collisions. Preliminary versions of this model have been 
presented in [io| . As it happens in most of the exist- 
ing codes, we shall work in the framework of the Dual 
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Parton Model (DPM) [ll| with string formation (valence 
strings, in the present case) based on standard Parton 
Distribution Functions (PDFs) with QCD evolution and 
string fragmentation via the Schwinger mechanism. The 
basic ingredients of the model are: 



• Formation of extended color fields or strings, mak- 
ing use of PDFs for valence quarks; 



Evolution with momentum scale Q 2 , as a conse- 
quence of the QCD evolution of the PDFs; 



• Fragmentation of each string with formation of a 
fast baryon (net-baryon) and other particles. 

With this simple model we show that a significant part of 
the physics of net-baryon production can be understood 
in terms of valence strings and Q 2 evolution, while sea 
strings play a less relevant role. 

When going from p-p to A-A collisions, it must be 
noted that the net-baryon results depend directly on the 
number of participants in the collision and thus on the 
collision centrality. The relation between the number of 
participants and the impact parameter of the collision 
can be established, for instance, in the context of the 
Glauber model [l^ . For A-A collisions, centrality thus 
plays an important role in our model. 

This paper is organized as follows. In section [Til the 
available experimental data and the predictions of cur- 
rent models are briefly reviewed and compared. In sec- 
tion [TTTJ our net-baryon model is described. The results 
of the model are presented and compared with data and 
with existing models in section [TV] A summary and brief 
conclusions are presented in section fVl 
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II. REVIEW OF CURRENT MODEL 
PREDICTIONS AND EXPERIMENTAL DATA 

The data presently available on net-proton or net- 
baryon rapidity is relatively scarce. The existing exper- 
imental results used in this study are briefly reviewed 
here. 

The most recent data are from the BRAHMS collabo- 
ration at RHIC [TH, [HI ; which published net-proton and 
net-baryon rapidity distribution for central Au-Au col- 
lisions at y/s = 200 GcV and yfs = 62.4 GeV. The 5% 
most central collisions were used at y/s = 200 GeV, while 
at y/s = 62.4 GeV the 10% most central collisions were 
selected. The procedure to obtain the net-baryon dis- 
tributions from the measured net-proton ones (properly 
accounting for strangeness and for neutrons) is detailed 
in the references. In the BRAHMS net-proton distri- 
butions, weak decay corrections (to remove the contri- 
bution from protons/antiprotons produced not at vertex 
but from strange baryon decays) were not applied. 

At the SPS (Pb-Pb collisions at y/s ~ 17 GeV), the 
NA49 collaboration published net-proton rapidity distri- 
butions with weak decay corrections included for cen- 
tral collisions (the 5% most central collisions were se- 
lected) [lj|. Net- proton results for different centrality 
bins have been presented by this experiment in (l6| . In 
this work data have been split into five centrality bins: 
0-5%, 5-14%, 14-23%, 23-31%, 31-48% and 48-100%. 

At lower centre-of-mass energies, net-proton rapidity 
distributions have been obtained at AGS (Au-Au at 
y/s ~ 5 GeV) [iTj for the most central collisions. Cen- 
trality cuts of 5%, 4% and 3% were used by the E917, 
E877 and E802 experiments, respectively. Weak decay 
corrections are negligible at AGS energies. 




tions and in terms of rapidity loss, defined as 

(Sy) =yp- (y) , (l) 

where y p is the beam rapidity and (y) is the mean net- 
baryon rapidity after the collision, given by 



(y) = 



N, 



2 f y " dn B ~ S (y) , 

y — i d v- 



part JQ 



dy 



(2) 



Here, N part is the number of participants in the collision 
and N B _g is the net-baryon number. 

Let us now turn to net-baryon production as imple- 
mented in the existing Monte Carlo models. QGSJET- 
II and EPOS Q are amongst the presently most 
widely used hadronic models in high energy and cosmic 
ray physics. To our knowledge, there is no systematic 
study comparing the predictions of these two models on 
net-baryon production between themselves or with ex- 
perimental data. 

In figure [TJ net-proton rapidity distributions obtained 
with QGSJET-II.03 and EPOS 1.61 for Pb-Pb collisions 
at y/s ~ 17 GeV are shown and compared with exper- 
imental data. According to [HI weak decay corrections 
have been applied to the data. Following [l8| , the impact 
parameter range (0 to 3.1 fm) corresponding to the cen- 
trality cut of 5% applied to data was selected in QGS JET 
and EPOS. While in EPOS an impact parameter cut 
is performed directly, in QGSJET this is implemented 
through a multiplicity cut Q- For EPOS, the results ob- 
tained leaving all particles free to decay and contribute to 
the net-baryon are shown together with those obtained 
switching off all weak decays. The 'all-decays' EPOS 
curve reproduces the trend seen in data, with an excess 
that could be due to the presence of net-baryon from 
weak decays. However, the strangeness effect compen- 
sation in the 'no-decays' curve seems to be too strong. 
For QGSJET-II, the 'no-decays' curve is shown and was 
found to be very similar to what is obtained from the 
model by default. It should be noted that in QGSJET 
the parton shower is fully developed before hadronisation 
is performed, and the connection to the experimental pic- 
ture is less direct. One should recall that QGSJET-II is 
not expected to perform very well at such low energies. 



FIG. 1: The net-proton rapidity distributions from EPOS 
1.61 (full and dotted lines) and QGSJET-II.03 (dashed line) 
for Pb-Pb collisions at y/s ~ 17 GeV are shown and compared 
with NA49 data (points) with weak decay corrections [l5| . 
For EPOS, the results obtained including (full line) and ex- 
cluding (dotted line) the strangeness contribution are shown. 
For QGSJET, the 'no-decays' curve is shown. 

The experimental results are usually presented in 
terms of net-proton and net-baryon rapidity y distribu- 



III. THE MODEL 

In the spirit of 0, H, [HI, [IH , the basic assumption of 
our model is that net-baryon production in proton-proton 
collisions is strongly correlated with the formation and 
fragmentation of two color singlet strings, each one with 
two valence quarks from one of the protons, and one va- 
lence quark from the other proton. This is schematically 
shown in figure 

Referring to the figure, let x\, a; 2 and X3 be the frac- 
tions of momentum carried by the valence quarks forming 
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FIG. 2: Schematic representation of a proton-proton collision, 
with the formation of two valence strings. 



string A. Quarks 1 and 2 are from the proton with pos- 
itive momentum in the proton-proton reference frame, 
and quark 3 is from the other proton. No transverse 
momentum is considered within this model. By choice, 
xi > x 2 > Xs, with < 0. The energy and momentum 
of each string are obtained adding directly the energy 
and momentum carried by each of the valence quarks. 
For string A: 



Estring = (%1 + X 2 + (~ aia))^ 

\/~S 

Pstring = (xi + X 2 — (—X3)) — 



y/{xl +x 2 ) (-x 3 ) s. 



string 



(3) 

(4) 
(5) 



where the quark momentum fractions X\, x 2 and x% are 
determined from the valence quark PDFs at an effective 
momentum scale Q 2 . For each y/s, an effective Q 2 de- 
rived from a fit to experimental data will be chosen (see 
section HV]) . In this work the CTEQ6M parton distribu- 
tion functions ppj were used. 

The simplest possible model for fragmentation is as- 
sumed. Each string decays into a baryon and a meson in 
the following way: the string is cut in two pieces and a qq 
pair is formed, from the vacuum, either between quarks 2 
and 3 (belonging to different protons) or between quarks 

1 and 2 (belonging to the same proton, with positive mo- 
mentum in the case of string A). The quark pair with the 
largest momentum difference is chosen. The string piece 
that inherits two valence quarks originates the baryon, 
whereas the string piece that inherits one valence quark 
originates the meson. This mechanism corresponds to 
the diagrams represented in figure [3] Diagram 1 corre- 
sponds to the case x\ — x 2 < x 2 — (—£3), in which quarks 

2 and 3 are chosen. In Diagram 2, x\ — x 2 > x 2 — (— X3) 
and string fragmentation occurs between quarks 1 and 
2. The weights of the two diagrams are, in this model, 
given only by kinematics. For string A the first diagram 
will be more probable (especially for low y/s). However, 
the weight of the second diagram can be as large as 40% 
above LHC energies. 

The qq pair formed from the vacuum was taken to be 
cither a mm or a dd, and the full quark combinatorics 
was then performed in order to determine the nature of 
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x\ - x 2 < x 2 - (-X3) 
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Diagram 2 

x\ - x 2 > x 2 - (-£3) 

FIG. 3: The two main valence string fragmentation diagrams. 
In the simple model assumed, the string is cut in two pieces 
and a qq pair is formed, from the vacuum, between the two 
quarks with the largest momentum difference. Diagram 1 
corresponds to the case X\—X2 < X2 — (—xs), in which quarks 2 
and 3 (belonging to different protons) are chosen. In Diagram 
2, x\ — X2 > X2 — (— X3) and string fragmentation occurs 
between quarks 1 and 2 (belonging to the same proton). 



the possible outcoming baryon. Both fundamental and 
excited states were considered, taking spin-dependent 
weights (2j+l). The decays of the unstable baryons were 
then performed and the outcoming nucleons included in 
the net-baryon calculations. The contribution from s 
quarks was not considered. The strangeness effect was 
estimated from the Schwinger model to be about 25% 
to 35%. The implications of this approximation arc dis- 
cussed below. 

The string mass and momentum distributions as a 
function of Q 2 are given in figure 0J The mass distribu- 
tion applies to both valence strings, while the momentum 
distribution for the case of string B is simply symmetric. 
The net-baryon rapidity distributions for different val- 
ues of Q 2 are presented in figure [SJ for proton-proton 
collisions at two different centre-of-mass energies. This 
figure illustrates not only the evolution with Q 2 but also 
the clear effect of kinematics, due to the increase of y/s, 
on the main features of the distribution. It should be 
noted that, as we are dealing with valence strings, when 
the Q 2 increases the mean rapidity decreases, simply due 
to PDF evolution. 

It is worth stressing that the inclusion of diagram 2, in 
addition to diagram 1, with weights determined simply 
by kinematics, reproduces some of the effects predicted 
in models with string junctions [U, [22j or popcorn [23[ 
mechanisms for the transport of baryon number from the 
beam rapidity into the central region y ~ 0. These effects 
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FIG. 4: Evolution with Q 2 (in (GeV/c) 2 ) of (a) the string 
mass and (b) the string momentum distributions for proton- 
proton collisions y/s — 200 GeV. The mass distribution ap- 
plies to both valence strings, A and B (see figure [2} • The 
momentum distribution is for the case of string A, and is 
symmetric for string B. 



can thus be achieved in a simple DPM model, based on 
valence strings and a Q 2 parameterisation. 

It is thus apparent that in this model the net-baryon 
rapidity dn/dy B ~ B in proton-proton collisions is ob- 
tained from the rapidity of the two baryons produced 
in the fragmentation of the valence strings, having the 
effective momentum scale Q 2 at each y/s as the only free 
parameter. 

Let us now consider the case of A-A collisions. As 
noted above, the net-baryon results depend directly on 
the collision centrality. The number of participants will 
thus play an important role in the model. We shall as- 
sume that the net-baryon rapidity in A-A collisions at a 
given y/s can be obtained from the net-baryon rapidity 
in p-p collisions at the same y/s computed at an effective 
Q\ value which depends both on y/s and on A, with a 
normalisation factor which is the number of participants 
per nucleus, N part /2: 

A B-B it 7 B—B 

f (Qi) * £!f± ■ p (Q 2 A ) , (6) 
ay a-a 2 ay p-p 

where both rapidity distributions are evaluated at Q\. 
As stated above, the model does not include ss pairs. We 



FIG. 5: Evolution of the net-baryon rapidity with Q 2 (in 
(GeV/c) 2 ) for proton-proton collisions at (a) y/s = 17 GeV 
and (b) yfs = 200 GeV. 



are thus assuming that strangeness does not considerably 
distort the distribution. 

When comparing to net-proton results, we consider 
that net-proton is roughly 1/2 of (B-B) [U]. As 
strangeness effects are not included in the model, the 
obtained rapidity distribution is well suited for compar- 
ing to data with weak decay corrections included. The 
strangeness effect is in this case simply a global factor e s , 
estimated from the Schwinger model to be about 25% to 
35%, which can be included in the normalisation factor: 

dn p ~ p N t dn B ~ B 

f (Q A ) *^=l.(l_ e .).£ (Q\) . 
ay A-A 4 ay p-p 

(7) 

Nuclear effects correction factors for the valen ce q uark 
PDFs were estimated using EKS98 [HJ and nDS [25| and 
found to be below 10-15%. These corrections are taken 
into account in the calculations. 

To address the energy evolution of the model, a re- 
lation between the effective Q 2 and y/s needs to be es- 
tablished. The effective Q 2 corresponds to the typical 
transverse size (area) of the parton (here, the valence 
quark). It is reasonable to assume, as in Regge phe- 
nomenology (2fjj . that the average number of partons in 
a nuclcon increases as a power of the centre of mass en- 
ergy y/s. Thus, R 2 h Q 2 ~ y/s/y/so, where Rh is the nu- 
cleon radius which we take as fixed. It then follows that 
Q 2 should grow, allowing for deviations from this naive 
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expectation, according to 

Q 2 = Q 2 o(0 "[GeV 2 ]. (8) 

Further, the effective A-A Q\ should be analogously re- 
lated to y/s, but account for an extra dependence on the 
number of participants. Since the number of partons in- 
volved in an A-A collision grows proportionally to the 
number of nucleons from each nucleus involved in the 
collision (i.e. oc N part /2) and the transverse size (area) 
of the nucleus grows proportionally to (N part /2) 2 / 3 , we 
can expect the typical inverse size of a parton Q 2 A to 
grow proportionally to . Also, we can ex- 

pect this simple geometrical estimate to be modified by 
nuclear effects and thus write for the A-A effective Q A , 
as a function of the p-p effective Q 2 : 

Q\ = (^f) a Q 2 - 0) 

The value of a will be estimated using the data at differ- 
ent centralities available for y/s ~ 17 GeV. 

Combining equations ([5]) and ©, we can thus write, 
for A-A collisions: 

Q\ = Ql(^f) C '[£)\GeV 2 ], (10) 

where the exponent A„ will be estimated performing a fit 
to the available data points. Below, the free parameters 
of the model will be fixed by adjusting the results of the 
model to the experimental data at each yfs and centrality. 

IV. RESULTS 

In this section, the available data on net-baryon and 
net-proton rapidity reviewed in section [ill will be used to 
fix the free parameters of the model. The ability of the 
model to describe these data in a consistent way will be 
evaluated and discussed. While the data at different cen- 
tralities allows to study the dependence on the number 
of participants and to address the relation between p-p 
and A-A, data at different centre-of-mass energies allow 
to establish and check the evolution with y/s of the model 
predictions. 

The experimental results are presented in terms of net- 
baryon rapidity in the case of BRAHMS and in terms 
of net-proton rapidity in the case of NA49 and of the 
AGS experiments. The model calculations arc performed 
using equation ((6|) for the net-baryon case and ([7]) for the 
net-proton one. 

The NA49 net-proton rapidity distributions for differ- 
ent centrality ranges at yfs ~ 17 GeV presented in [lrj | 
were used in the first step of this analysis. Our model was 
fitted to these data points in each centrality bin, taking 
as free parameters the effective Q 2 A and a normalisation 
factor n. The data points and the results of the model are 
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FIG. 6: The results of the present model for net-proton rapid- 
ity are compared to NA49 data in different centrality ranges 
(from top to bottom: 0-5%, 5-14%, 14-23%, 23-31%, 31-48% 
and 48-100%.) [H. 



shown in figure [5J A satisfactory agreement is achieved 
for all the centrality ranges. The apparent difference be- 
tween data and the model in the mid-rapidity region, 
could be due to the fact that the experimental error bars 
do not include the error associated to weak corrections. 
On the other hand, difractivc effects arc not included in 
the model. This accounts for the poorer description of 
the high rapidity extremes in peripherial collisions, where 
such effects are expected to become more relevant. It is 
worth noting that the effect seen in data displacing the 
maxima of the distribution to higher rapidity values as 
the centrality of the collision decreases, is reproduced in 
the model through a decrease of the fitted Q 2 A value. 




N pa „ 



FIG. 7: Results of the fit to the NA49 rapidity distributions 
at different centralities: the fitted and expected value of the 
normalisation factor n are plotted against the value of N par t 
computed in [l(| for each centrality bin. 



The normalisation factor n obtained from the fit can 
be compared with what would be expected from equa- 
tion |(7J) using the number of participants estimated by 
the experiment in (lH | and assuming e s = 0.25. This 
comparison is shown in figure where the fitted and ex- 
pected values of n are plotted against the value of N part 
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computed in the reference for each centrality bin. We see 
that a good agreement is found. The relation between n 
and Np art is thus well understood, giving us confidence 
on the model and on the fitting procedure. 
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FIG. 8: Results of the fit to the NA49 dn/dy distributions 
at different centralities: the obtained (Q'a, N par t) points are 
shown and fitted with equation © (full) and with a straight 
line (dashed). 

We thus have, for each centrality, a number of partic- 
ipants and a value of Q\. These results can be used to 
estimate the exponent a on equation ©, which relates 
the effective Q 2 in p-p and A-A collisions. This was done 
by fitting equation ([9]) to these {Q\,N part ) points, leav- 
ing Q 2 and a as free parameters. The data points are 
shown in figure [HI together with the fit result, which cor- 
responds to a = 0.53^0 i3- A linear fit to these points 
is also plotted, showing that the limited range covered 
in N part does not allow to strongly constrain the form of 
the function describing this dependence. The constant 
term added to Q 2 at N part =2 in the case of the linear 
fit is however hard to motivate. It should be noted that 
a simple interpolation between the points is sufficient to 
relate Q 2 A and Q 2 in the range of N par t relevant for the 
interpretation of the existing data. The exact form of the 
function would however become relevant when extrapo- 
lating to p-p collisions. The uncertainty on the relation 
between Q 2 and Q 2 A is thus at present too large to allow 
such an extrapolation. 

Turning now to the evolution with y/s, we consider the 
results available for central collisions at RHIC, the SPS 
and AGS QJ, Q GJ, M, E3 ■ As above > in order to fully 
define the model we need to choose an effective Q 2 A value 
and a normalisation factor related to N part (and to e s , 
in the case of net-proton). In practice, this was done by 
fitting the model predictions to the existing rapidity data 
using as free parameters the effective Q 2 (and assuming 
the relation to Q 2 A we have just inferred from centrality 
data) and a normalisation factor n. The data points and 
the results of the fit are shown in figure [HJ As stated 
above, the experimental results arc presented in terms of 
net-baryon rapidity in the case of BRAHMS and as net- 
proton rapidity in the case of NA49 and of the AGS ex- 
periments. The model calculations were performed with 
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FIG. 9: The results of the present model for: (a) net-proton 
and (b) net-baryon rapidity are compared to experimental 
data from central A-A collisions at different centre-of-mass 
energies. See [H, [f4|, EH, E3 f° r details on the data points. 



equation (O for the net-baryon case and ([7]) for the net- 
proton one. A good description of the existing data is 
achieved. At y/s ~ 5 GcV only the points up to the 
nominal beam rapidity were considered, since our model 
does not include low energy effects relevant at these ener- 
gies, such as Fermi momentum effects, and therefore has 
no mechanism to reproduce these data. At all the other 
centre-of-mass energies all data points were included in 
the fit. The present RHIC data at y/s = 200 GeV cover 
only the mid-rapidity range, leaving the fit largely uncon- 
strained. For this reason, at this centre-of-mass energy 
the normalisation factor was fixed to the value estimated 
in reference [l3[ and only Q 2 was left as a free parameter. 
In all the other cases the normalisation was left free. 

The values obtained for the effective Q 2 , the normalisa- 
tion factor n and the number of participants N part (com- 
puted from equations ((6|) and (O) are given in table [H 
where the number of participants estimated by the exper- 
iments [III, 0, [H, QjJ f° r the relevant centrality region is 
also given. The fitted N par t values are in good agreement 
with the expectations. The discrepancy at y/s = 5 GeV 
is due to the effect in the normalisation factor of exclud- 
ing the points above the nominal centre-of-mass energy 
from the fit. 

As discussed above, strangeness contribution should be 
simply a global factor when comparing with net-proton 
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FIG. 10: The effective Qa values at different centre-of-mass 
energies chosen by tuning the model to the experimental data 
are shown. The line shows the fit using eq. (|10[) and consid- 
ering Au-Au collisions. The shaded areas correspond to la 
variations of the fit parameters. The inset figure shows the 
extrapolation of the fit to higher energies, covering the LHC 
and very high energy cosmic rays. 
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Collision 
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66.6j£* 


266.4l? 2 6 4 


344 ±6 


17 


Pb-Pb 


7R+ 05 
u - '°-0.07 


67.7111 


361.1±?;| 


362 ± 12 
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Au-Au 


o.77tg:g 


148.21™ 


296.4±gS:g 


314 ±8 


200 


Au-Au 


1 14+0.29 
- 1 — l<l -0.20 






357 ±8 



TABLE I: Results of the fit to the effective Q 2 and the normal- 
isation factor n at the different centre-of-mass energies. The 
number of participants obtained from equations ((6]) and J7} 
(N P art) and estimated by the experiments (N^/ t ) are also 
given. 




FIG. 11: Predictions of the present model for (a) net-baryon 
rapidity for central Pb-Pb collisions at y/s = 5.5 TeV and (b) 
rapidity loss as a function of y/s. The full line is the prediction 
of the model for Au-Au collisions with a 5% centrality cut. 
The experimental results available at the different centre-of- 
mass energies are also shown (see key inside the figure), and 
compared to the predictions of the model for the same type 
of collisions and centrality cut (stars). 



results and was taken to amount to e s = 0.25 at SPS 
energies and to be negligible at AGS. It could slightly 
distort the spectrum in the comparison with net-baryon 
rapidity distributions at RHIC, but this is not yet rele- 
vant given the error bars of the available data points. 

Assuming equation (TIT)]) to describe the evolution of 
the effective Q\ with centrality and i/i, the exponent A„ 
was determined by fitting this equation to the (s/s, Q 2 ) 
points in table [U This is shown in figure 1101 where the 
points are the Q 2 values adjusted above, the line is the 
fit to these points with equation (flO|) and the shaded 
areas correspond to la variations of the fit parameters. 
The value obtained for the exponent was A„ = 0.25!toiij 
taking y/s^ ~ 17 GeV. 

The predictions of the present model for net-baryon 
rapidity and rapidity loss at higher centre-of-mass ener- 
gies were then obtained and are shown in figure [TT1 The 
width of the curves corresponds to varying the Q 2 within 
the la band shown in figure ITTJl In figure [TTT a) the net- 
baryon rapidity distribution in central Pb-Pb collisions 
at y/s = 5.5 TcV is presented. Following reference [27| . 
the number of participants corresponding to selecting the 
6% most central collisions was set. In figure HUb) the 



available data points on the rapidity loss are shown and 
compared to the predictions of the model. For the model, 
two types of results are shown: the full line is the pre- 
diction for Au-Au collisions with a centrality cut of 5%, 
while the stars show the prediction directly comparable 
with the experimental data point at the same y/s (same 
type of collisions and centrality cut). A good agreement 
is found in all cases. The fact that the la band is quite 
narrow in this variable when compared to the Q 2 one re- 
sults from the weak dependence on Q 2 reached at high 
Q 2 values (and high energies). 

Finally, the fraction of the centre-of-mass energy car- 
ried by the net-baryon as a function of y/s was computed 
and is shown in figure [12] (together with the la bounds) . 
Theprediction is for Au-Au central collisions. According 
to 0, RHIC data indicate about 27% of the initial en- 
ergy remaining in the net-baryon after the collision. This 
result is also shown in the figure. 

At RHIC energies, the present model is somewhat be- 
low the measured value, but within 2a. At higher en- 
ergies, a sizable amount of energy is still associated to 
the net-baryon. It should be noted that high energy ef- 
fects such as string percolation may change these prcdic- 
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FIG. 12: Evolution of the fraction of energy carried by the 
net-baryon with y/s. The prediction of the present model 
(with a shaded band corresponding to the la variation of the 
fit parameters) for Au-Au collisions with a 5% centrality cut 
are shown and compared to the RHIC estimate given in [fj. 



tions [2. 



V. SUMMARY AND CONCLUSIONS 

A simple and consistent model for net-baryon produc- 
tion in high energy proton-proton and nucleus-nucleus 
collisions was presented. The basic ingredients of the 
model are valence string formation based on standard 
PDFs with QCD evolution and string fragmentation via 
the Schwinger mechanism. The results of the model were 
presented and compared with data from RHIC, SPS and 
AGS. 

The free parameters in the model, the effective Q 2 and 
the number of participating nucleons, were fitted to net- 



proton data published by AGS experiments and by NA49 
and at yfs ~ 5 GeV and -y/s ~ 17 GeV, respectively, and 
by the BRAHMS experiment and ^/s = 62.4 GeV and 
= 200 GeV. A good fit to the data is obtained within 
this model, with the number of participants matching 
the calculations in the literature. The results show that 
a good description of the main features of net-baryon 
data is achieved on the basis of this simple model, in 
which the fundamental production mechanisms appear 
in a transparent way. 

Using the available net-baryon data at different ccntre- 
of-mass energies and centralitics, a relation between the 
effective momentum scale and the centre-of-mass energy 
was motivated and a prediction was obtained and extrap- 
olated to higher energies for the evolution with yfs and 
Npart of the fraction of the initial energy carried away 
by the net-baryon. A sizable amount of energy may be 
associated to the net-baryon, even at high energies. 
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